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ABSTRACT 

The existence of large scale voids in several galaxy surveys suggests the occurrence 
of an inflationary first order phase transition. This process generates primordial bub- 
bles that, before evolving into the present voids, leave at decoupling a non-Gaussian 
imprint on the CMB. 

In this paper we evaluate an analytical expression of the collapsed three point 
correlation function from the bubble temperature fluctuations. Comparing the results 
with COBE-DMR measures, we obtain upper limits on the allowed non-Gaussianity 
and hence on the bubble parameters. 
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1 INTRODUCTION 

In the recent past the number of papers devoted to non- 
Gaussian anisotropies on the Cosmic Microwave Background 
(CMB) has increased dramatically. This new investigation 
field is, in fact, a powerful tool to distinguish between the 
theories of structure formation based on inflation and those 
based on topological defects. Quantum fluctuations pro- 
duced in inflationary models are scale invariant and have 
a Gaussian distribution. Thus we expect that three-point 
correlation function of the CMB temperature vanish (Falk 
et al. 1993; Luo & Schramm 1993; Gangui & al. 1994). 

On the contrary in models with topological defects the 
primordial density perturbations are scale dependent and 
non-Gaussian (Avelino et al. 1998): hence we expect some 
deviations from Gaussianity in higher order correlation func- 
tions. In this context we may also include the extended infla- 
tion model (La & Steinhardt 1989), because during the in- 
flationary epoch we have a first order phase transition, that 
generates bubbles of true vacuum. These voids contribute to- 
gether with ordinary quantum fluctuations to structure for- 
mation. This possibility has been investigated (Occhionero 
& Amendola 1994; Amendola et al. 1996): it has been shown 
(Occhionero et al. 1984, 1997) that primordial bubbles may 
be associated with the observation of large scale voids in sev- 
eral galaxy surveys (Kirshner et al. 1981; de Lapparent et al. 
1989; da Costa et al. 1996; El Ad, Piran & da Costa 1996, 
1997). Since these defects can also produce non-Gaussian 
anisotropies on the CMB, we may obtain some limits on 



the bubble parameters comparing observations with non- 
Gaussian predictions. 

So far differents statistical tests have been applied to 
COBE-DMR sky maps (Kogut et al. 1996) and the results 
have been in agreement with the Gaussian models. Although 
recently two groups have detected non- Gaussian signal in 
COBE data (Ferreira et al. 1998; Pando et al. 1998), sub- 
sequently this has been shown to derive from a systematic 
effect in the data (Banday et al. 1999). On the other hand, 
Bromley & Tegmark (1999) tried to argue that the COBE 
4 year data were in fact Gaussian. 

In this paper we compare the level of non-Gaussianity 
produced in bubble models with the COBE data. We evalu- 
ate the three point correlation function in CDM models that 
contain also primordial bubbles. Comparing the numerical 
results with the COBE-DMR measures (Kogut et al. 1996; 
see also Hinshaw et al. 1994, 1995) we obtain upper limits 
on the parameters of the voids in agreement with galaxy 
surveys observations. 

2 METHOD 

The imprints of bubbles on the CMB has been studied in 
several papers (Baccigalupi, Amendola & Occhionero 1997 
; Amendola, Baccigalupi & Occhionero 1998; Baccigalupi & 
Perrotta 1999). The presence of the primordial voids induces 
a Sachs- Wolfe effect and an acoustic perturbation propagat- 
ing up to the sound horizon on the photon distribution. As 
a consequence, the induced temperature fluctuations of in- 
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dividual bubbles are composed of a central spot and some 
concentric hotter isothermal rings; this pattern has been cal- 
culated by numerical integration of the Boltzmann equations 
in Amendola, Baccigalupi & Occhionero (1998). It is found 
that the bubble signal depends on the radius R of the void 
and on the central (negative) density contrast S. We shall 
distribute N voids on the CMB sky and use the fraction X of 
the space that the voids fill today as a free parameter, where 
X = NR :i /3L'iAL h (Amendola et al. 1998) where L h is the 
horizon radius and ALh is the thickness of the last scat- 
tering surface. We consider bubbles of size R = 30/i^ 1 Mpc 
at decoupling: due to their overcoming growth, these voids 
have today radii around 20 ~ 60/jT 1 Mpc, like those ob- 
served in galaxy surveys (da Costa et al. 1996; El Ad et al. 
1996; 1997). 

In simulated COBE maps, due to the low resolution of 
the satellite, the signal of the individual bubbles looks like 
dark spots confused amidst Gaussian anisotropies; their ef- 
fect appears only when caculating the correlation functions 
of maps containing many bubbles. The temperature fluctu- 
ation may be decoupled in two terms: 



A {0, <p) = Ac (e,<p) + A v (e,v)- 



(i) 



The first term Ago» S s (0, ip) is the Gaussian temperature 
fluctuation field produced by the primary anisotropies; the 
second term is the voids signal, that vanishes in directions 
where there are not bubbles. In order to compare the pre- 
dictions of the model with experimental data, we calculate 
the collapsed three-point function 
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where W l is the window function of the experiment, Pi are 
the Legendre polynomials, a™ are the multipole coefficients 
of the spherical harmonic expansion and where 
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with ^ ^ ) the Wigner 3 J symbol. We compare Ci{a) 
to the pseudo three-point collapsed function of Kogut et al. 
(1996): obviously the two are identical due to the absence of 
noise in our case. Since the Gaussian term and the signal of 
the bubbles are not correlated, we may write the three-point 
correlation function as sum of two separate contributions: 



r^Gauss . 
3 = 3 + O3 . 



(4) 



The contribution to C3 from gaussian fluctuations is not 
zero. This contribution may arise from non-linearities in 
the inflationary dynamics or from non-linear growth of the 
perturbations. However, using the analytical expression for 
the C§ auss (a) computed in Gangui et al. (1994) it can 
be seen that the level of non-Gaussianity produced by the 



non-linearities in the inflation dynamics is smaller than 
that arising from the non-linear growth (Mollerach et al. 
1995), and that the latter is much smaller that produced 
by the bubbles. In fact, on the angular scales probed by 
COBE-DMR, Mollerach et al. (1995) found an amplitude 
(C 3 H ~ s (a)) ~ 0.1 fiK 3 , while we find that the contribu- 
tion of the voids is larger by several orders of magnitude: 
(C 3 v (q)) ~ 10 4 fiK 3 . Therefore, we neglect C§ auss in the 
following. 

To calculate C3 (a) we use the same approch of texture- 
spot anisotropies (Magueijo 1995, Gangui & Mollerach 1996, 
1997). The temperature fluctuations produced by a random 
distribution of bubbles, in the 7 direction, is simply the 
superposition of the signal of all the bubbles, and can be 
written as Ay (7) = ~}2 n b n f n (-y n ,a). Here, the signal of 
the n-th bubble has been decomposed as a overall ampli- 
tude b n (corresponding to the central temperature fluctua- 
tion) and a density profile /n(7 n ,«) where a is the angle 
measured from the bubble center. The dependence on the 
parameters R n and 8 n is contained only in the amplitude; 
it is to be expected that this dependence is linear in S n 
and quadratic in R n , since the central temperature fluctua- 
tion is dominated by the Sachs- Wolfe effect. The expression 
b n = 8 n (R n /20H~ 1 ) 2 is indeed an accurate fit in the range 
we are interested (see Amendola et al. 1998). The profile 
fn{j n ,a) contains the full effect of the acoustic oscillations 
and the adiabatic fluctuations, and has been obtained nu- 
merically in Baccigalupi & Perrotta (1999). 

We expand Ay in spherical harmonics and obtain the 
multipole coefficients 

°r = 2r^E^'(u ( 5 ) 

n 

where F l n is the Legendre trasform of the intensity profile, 
F n = J dn a /„( 7n ,a)P ( (cosa). (6) 

Inserting (5) in (2) the collapsed function reduces to: 
Cf(a) = 4tt J2 Ph(cosa)W ll W h W l3 J hl213 

x 
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We take the window function of COBE to be ~ e 



-l(l + l)a 2 /2 



with a = 3.2°. We assume now that there are N identical 
voids on the CMB sky. Developing the sum on m, ni and 
ri3 we obtain three terms that represent the contribution to 
the C3 (a), when the bubble signals are not correlated and 
when are correlated two by two or three by three and etc. 
We take into account the correlation at lowest order, in other 
words we consider just the first two terms. Then the mean 
value of (7) for a Poissonian bubble distribution on the sky is 
obtained substituting the sum on the bubble index with an 
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integral over the whole sky. In fact, the number of bubbles 
in a circular ring centered on a single bubble is proportional 
to angular extension of the ring, therefore we have: 
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F h {9)= J f{0 + a)Pi 3 (cos a)d(cos a). (11) 

Using the same approch, after a tedious calculation, we 
have found an analitycal expression for the variance a\[a) = 



that we do not report for shortness. 



We compare the experimental data with the behaviour of 
the (C$ (a)) for differents values of the parameters 8 and 
X. When (C^(a)) ±ov(a) is larger than COBE data plus 
the noise and cosmic variance, we have some constraints on 
the parameters of our model. 



3 RESULTS 

The COBE data has been taken from Kogut et al. (1996). 
We assume a fraction of bubbles corresponding to 0.31 < 
X < 0.54, consistent with da Costa et al. (1997). We have 
computed the C^(a) for 0.001 < 8 < 0.0026, without 
dipole and quadrupole contribution, / m i n = 4. In the fig- 
ures we report the behaviour of the C% (a) for two values 
of 8 = 0.002, 0.0012. The oscillating behaviour of the plots 
is due to the sum of the Legendre polynomials in (9). In 
the plots the errorbars are the a(a)'s. The level of the cos- 
mic variance a (a) generated from the model is very high 
for a < 40°, while it is small on the large angular scales, 
a > 45°, where the contribution of the lowest multipoles is 
small. In figure (1) we have the model with 8 — 0.002: we 
may note that for X — 0.54 the signal is larger than cosmic 
variance and the observed data points, while X — 0.31, the 
plot is marginally consistent with the experimental data. In 
figure (2) we report the (C¥(a))T$ for 8 = 0.0012: it fits 
the COBE data very well. Notice that in the range a > 50° 
the behaviour of the collapsed function seems to follow the 
trend of the COBE measures. Values of 8 < 0.0012 produce a 
C$ (a) within the cosmic variance band and smaller than the 
COBE data. In this case the observations do not impose con- 
straints and we may obtain only an upper limit on the value 
of 8. We have applied a \ 2 analysis to our models. In figure 
(3) we report the confidence regions with a confidence level 
set to 99.9% (grey region) and to 99.5% (black region). We 
may note that all models with 8 > 0.0017 are ruled out by 
the experimental data. Then we may conclude that although 
the bubbles produce a non-Gaussian signal on the CMB, this 
is in agreement with the present observation provided that 
the density contrast 8 < 0.0017 or X < 54%. So we obtain 



Figure 1. The points arc the COBE data while the thick lines 
are the cosmic variance of a Gaussian random field (Kogut et al., 
1996). The plots are models with <5 = 0.002 and X = 0.54 (dashed 
line) and X = 0.31 (thin line). The errorbars are the variance of 
our models. 



Figure 2. Same as in Fig. 1, but now <5 : 

observations is now obtained. 



0.0012: agreement with 



a constraint stronger than that found in Amendola et al. 
(1998), where the bubble power spectrum was compared to 
the measures of the CAT experiment. The next high reso- 
lution experiments, like MAP and Planck, and the recent 
observations of Boomerang and Maxima should be able to 
detect the voids signal on the CMB. In fact these missions 
can probe the multipoles I > 100, where the contribution of 
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Figure 3. Confidence regions for different values of <5 and X, the 
confidence level is set to 99.9% (grey region) and to 99.5% (black 
region). Models with 5 > 0.0017 are ruled out by the COBE 
measures. 



the bubbles is important, and the effects on Cs(a) may be 
large. 



4 CONCLUSION 

Several galaxy surveys found huge spherical voids in the 
matter distribution, the galaxies lying in the surrounding 
shells: these structures may be generated in inflationary 
models with first order phase transitions. These bubbles 
produce a non- Gaussian signal on the CMB. We analyse 
this signal developing an analytical expression for the three- 
point collapsed function of a bubble distribution, using the 
formalism of Magueijo (1995). Our free parameters are the 
density contrast and volume fraction of the bubbles, while 
the radius R is fixed to a value consistent with the galaxy 
surveys. We compare the behaviour of the three point col- 
lapsed function for the bubble model with the COBE data. 
We obtain a constraint on the value of 8: in fact, the ex- 
istence of the voids at decoupling is not in contrast with 
the measures of the COBE three-point collapsed function, 
provided 6 < 0.0017 or X < 0.54. This still leaves plenty 
of room for the bubbles to cooperate efficiently to structure 
formation, both via the central voids and via the possibil- 
ity of shocking on the outher shell: in fact a central density 
contrast of 0.001 can still evolve linearly in an empty void 
by today. More information will be obtained comparing the 
results of the future high resolution experiments. 
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